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ABSTRACT 

We derive from first principles equations governing (a) the quadrupole tensor 
of a star distorted by both rotation and the presence of a companion in a possibly 
eccentric orbit, (b) a functional form for the dissipative force of tidal friction, based 
on the concept that the rate of energy loss from a time-dependent tide should be a 
positive-definite function of the rate of change of the quadrupole tensor as seen in 
the frame which rotates with the star, and (c) the equations governing the rates of 
change of the magnitude and direction of the stellar rotation, and the orbital period 
and eccentricity, based on the concept of the Laplace-Runge-Lenz vector. Our analysis 
leads relatively simply to a closed set of equations, valid for arbitrary inclination of 
the stellar spin to the orbit. The results are equivalent to classical results based on the 
rather less clear principle that the tidal bulge lags behind the line of centres by some 
time determined by the rate of dissipation; our analysis gives the effective lag time as 
a function of the dissipation rate and the quadrupole moment. 

We discuss briefly some possible applications of the formulation. 

Subject headings: stars: binary, triple 

1. Introduction 

We present a partly novel derivation of some basic results, and some new results, for the 
'equilibrium tide' model of tidal friction (Hut 1981). We derive the force and couple on a binary 
orbit that arises from dissipation of the equilibrium tide, starting only from the principles that (a) 
the rate of dissipation of energy should be a positive definite function of the rate of change of the 
tide, as viewed in a frame which rotates with the star, and (b) the total angular momentum is 
conserved. This model leads to a force and couple which are similar to those derived classically on 
the assumption that the tidal bulge lags the line of centres by some small time that is determined 
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by the timescale of dissipation. However in our model the lag time is determined to be proportional 
to the quadrupole moment, rather than assumed to be independent of it. There is a clear relation 
between the rate of dissipation and the lag time. We also determine the relation between the 
hypothesised dissipation constant and the turbulent viscosity within the star. However there 
remains the more difficult task of determining the turbulent velocity field. 

There are two principle formalisms for discussing the effect of tidal dissipation on an 
orbit: the 'equilibrium tide' model (Darwin 1880, Alexander 1973, Hut 1981, Hut 1992) and 
the 'dynamical tide' model (Fabian, Pringle & Rees 1975, Zahn 1975, Press Sz Teukolsky 
1977, McMillan, McDermott & Taam 1987, Goldreich & Nicholson 1989, Lai 1996, Kumar & 
Goodman 1997, Kumar & Quataert 1997, Papaloizou Sz Savonije 1997, Savonije et al. 1983, 
1984, 1995, 1997, Mardling 1995). In the former the star is perceived as being perturbed to a 
shape closely approximating the shape it would have in equilibrium with the (time-dependent) 
gravitational-centrifugal potential of the system. Dissipation by viscosity of the internal motion 
caused by the time- variation of the tide causes the orbit to tend towards circularity, and causes 
the spin to be brought into parallelism with the orbit, until the tide ceases to be time-dependent 
in the limit that the orbit is circular and the spin parallel to the orbit (Hut 1980). In the latter 
formalism, the star is perceived as an oscillator with a number of normal modes that are excited 
by the time-dependent potential as the companion passes by at periastron. In a close encounter 
from an initially parabolic (or highly eccentric) orbit, energy is extracted from the orbital motion 
and locked up in the oscillation of the star's normal modes. Each mode is damped at some rate, 
the higher-order modes being damped more rapidly. If the modes are all damped during the long 
apastron phase, then the process can repeat itself at the next periastron, and thus represents a 
long-term drain of orbital energy intothermal energy. Since the angular momentum is constant 
while energy is dissipated, we expect the orbit to become more circular. If the modes damp 
slowly, they may still be active at the next periastron, and because the frequencies of oscillation 
are not necessarily commensurate with the frequency of the orbit, each periastron passsage may 
be different from the previous one; some oscillations may be in phase and excited resonantly, while 
others may be out of phase. It is possible for the orbit to be decircularised instead of circularised, 
and in some circumstances it may vary chaotically. Even for circular orbits, dynamical tides can 
continue to operate, as long as there still are deviations from synchronous rotation (Zahn 1975, 
Goldreich & Nicholson 1989). 

These two models are not necessarily in competition. Rather, it is likely that the first prevails 
in cases where an orbit is only moderately eccentric, and reasonably wide in the sense that both 
stars are small compared with their periastron separation. The latter may prevail in circumstances 
('tidal capture') where two stars approach initially on a parabolic or slightly hyperbolic orbit, 
and interact briefly and dramatically near periastron. In the former case it is likely that only the 
lowest order (quadrupolar) mode will be significantly excited, as in the case of tides on Earth. We 
concentrate here on the equilibrium-tide model, which we develop using a mathematical-physical 
treatment which we believe to be unusually clear and simple. 
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In § 2, we derive an approximation to the quadrupole tensor of the mass distribution in a 
star of mass rrt\ distorted by the combination of (a) uniform rotation J), and (b) the presence 
of a companion of mass at position d(t) relative to the centre of *1. SI is not necessarily 
perpendicular to d, and neither is d. We deal only with the lowest, i.e. quadrupole, mode because 
it is much the simplest; we can operate entirely within the linear regime, and extract results that 
are exact to first order. Including higher modes involves ultimately working in the non-linear 
regime, and one gains only a modest increase of accuracy, if any, at the expense of considerable 
extra complexity In § 3, we show that if the mechanical energy of the motion (orbital energy plus 
rotational energy) is dissipated at a rate which is proportional to the square of the rate of change 
of the quadrupole tensor, then a force (and associated couple) is introduced into the system which 
is equivalent to the assumption that the tidal bulge lags behind the line-of-centres by some time 
which depends on the rate of dissipation. In § 4, we use the concept of the Laplace-Runge-Lenz 
(LRL) vector to determine relatively easily the effect on the binary orbit of an arbitrary perturbing 
force. In § 5 we apply this analysis to the dissipative force obtained in § 3, in combination with the 
non-dissipative force that is already present simply because of the quadrupolar distortion of § 2. 
The standard equations are recovered relatively simply in the case that SI is parallel to the orbital 
angular momentum; and we obtain the complete equations for the case that they are not parallel. 
In § 6 we derive the tidal velocity field and its rate of dissipation by turbulent viscosity. This leads 
to a specific value for the dissipation constant of § 3. In § 7 we determine briefly the meridional 
circulation current (Eddington 1925, Vogt 1925, Sweet 1950, Mestel 1965) that is implicit in the 
distortion determined in § 2. However this circulation current must be taken with a pinch of salt 
since it shows several singularities in the outer layers of hot stars. In § 8 we discuss uncertainties 
in the model, and mention briefly possible applications of the model to (a) the circularisation of a 
pulsar orbit around a normal star, and (b) the effect on a close binary orbit of the presence of a 
distant third star in a hierarchical triple system. 



2. The Equilibrium Tide 

Consider a binary consisting of an extended star (*1), which rotates with constant uniform 
angular velocity SI, and a point-mass companion (*2) at a distance d(t) from it. In a frame which 
rotates with *1, and has origin at the CM of *1, the equation of motion of the fluid of *1 is 

Dv 1 , 1_ ,o Gni2 Gwi2 r.d . ,. 

— + 20Av = — Vp-V0 , 4> = fa- -SlAr 2 --—j- j- , (la, 6) 

Dt p 2 |r + d| a 6 

where r is a vector from the CM of the extended star to a general point in it, and D/Dt is the 
derivative following the motion. The potential <f> is the combination of the potential of self-gravity 
fa, where 

V 2 fa = AttGp , (2) 

along with centrifugal force, and with the tidal potential of the companion. The last term in Eq. 
(lb) allows for the fact that the frame centred at the CM of *1 is accelerated relative to an inertial 



-4 - 



frame; it cancels the first-order term in r/d from the second-last term. To lowest order in r/d, 

cf> « fa - hl 2 r 2 + ifiVP 2 (cos 9) - ^f^P 2 ( C os 9') , (3) 

where 9 is the polar angle measured from Q, and 9' is the polar angle measured from d. 

If the spin and orbit are not aligned, and/or the orbit is elliptical, the star cannot be entirely 
static, i.e. v ^ 0, but we assume that v is small so that at some level of approximation pressure 
balances rotational and gravitational forces: 

Vp = -pV<p . (4) 

Hence p, p are constant on equipotentials. We can then define variables V, m, r* - volume, mass 
and 'volume radius' respectively - which are also constant on equipotentials. V is the volume 
contained within an equipotential and m the mass, and they relate to r* by 

f r. s = V{4>) , ^ = txrM*) ■ (5) 

We can see that, since the distance between adjacent equipotentials along a normal is SI = S(j)/\V(f)\, 

dV [61 [ 

where dT, is an element of area. 

The constancy of p, p on equipotentials, and hence of T also at least in the case of chemically 
homogeneous stars, means that in the heat equation we have to introduce a velocity field, the 
'circulation current', to maintain thermal equilibrium (§ 7). We assume this field is sufficiently 
small that we can still neglect Dv/Dt, Jl A v in Eq. la. The theory developed in §§ 2 - 6 does 
not require any knowledge of the circulation velocity, only supposing that it is small so that the 
approximation of p, p constant on equipotentials is valid. 

The following subsections (i) - (viii) develop the following points as briefly as possible: 

(i) a definition of the quadrupolar-distortion parameter a(r) 

(ii) the quadrupole moment as an integral of a 

(iii) the surface boundary condition for a, in the case of rotation but no companion 

(iv) the second-order DE satisfied by a in the interior 

(v) analogous to (iii), but for companion and no rotation 

(vi) an alternative form of the surface boundary condition (Schwarzschild 1958) 

(vii) an illustrative but too crude approximation for a 
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(viii) the potential force fQD (in addition to point-mass gravity) between bodies when one has a 
quadrupolar distortion due to both rotation and a companion. 

(i) We consider only the lowest-order, quadrupolar, distortion. An equipotential surface r = r(r*, 6) 
will be approximated by 



r « r*{l - a(r*)P 2 (cos 9)} 



and conversely r* ps r{l + a(r)i-2(cos 9)} , (7a, b) 



where 9 is the angle from the axis of symmetry. Since a is first order, it can be thought of as a 
function of either r or r*. 

(ii) The contribution to the quadrupole moment q from mass between equipotential surfaces 4> and 
4> + dcj) is 



— = pr F2 (cos 1 



d£ 



2vrr 2 sin 0d0 



r 2 P 2 (cos6») , 



(8) 



where we use Eq. (6) to estimate the volume element. Now, is a function of r* only, and r* is a 
function of r, 9 given by Eq. (7b), so that 



r*) 



dr* 



dr„ 



dra(r) 

1 H -. -T2 

dr 



dr* 



drMn) 

1 H -. F-i 

dr* 



Hence 
dq 

dr* 

and so 



p(r* 



dr*a(r*) 
dr* 



r* 4 (l -4aP 2 )^2 2vr sin 9d9 



An 



Aa + 



dr*a 



dr* 



fnii 

/ (5a + r*a')r* 2 dm 
Jo 



1 f mi 

'5 



(9) 



(10) 



(11) 



(iii) The self-gravity term 0i of Eq. (3), outside the star, is just the sum of the monopole and 
quadrupole terms, and so the centrifugal-gravitational potential distribution outside the star is 



Gmi + Gq^ & + 1 v _ 1^2 r 2p 2(cos( 
r r J 3 3 



If ri is the radius of *1, this means that on the stellar surface (r* = ri), using Eq. (7a) 



Gmi 1_ 2 2 / Gmi Gg rot l o2 2\ n / m 
+ -ft 2 r 2 + ai + — \ -^ r{ P 2 (cos0) 



where ai = a(ri). This must be independent of 9, so that 



ai 



Using Eq. (11), this leads to 

n 2 rf 1 



ai 



jot 



n 2 r 5 



n 2 



,,TOt 



2>Gm\ m\r\ 
, Q 



(12) 



(13) 



(14) 



3Gm 1 1-Q 



3G 1-Q 



If r 2 dm(5a + ra') 

with Q = v 2 ' 

5 mir^a(ri) 



(15a, 6, c) 
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Once again, since a is first order we do not need to distinguish between a(r*) and a(r) in the 
integral. 

(iv) Evidently to determine Q we now have to determine the functional form of a(r). From Eq. 
(7b) throughout the interior, rather than just at the surface, we can obtain V 2 r* and |Vr*| 2 to 
first order in a and its derivatives, and since <f> as well as p is a function of r* only we find that 



V 2 6 = <h"\VrJ 2 + d)'V 2 r* = 6" + —6' + 



2a 

(t>'(ra" + 4a'--^) + 2<p"(ra' + a) 



b'Po 



= 4vrGp(r*) - 2Q 2 . (16) 

Outside the square brackets a prime means a derivative w.r.t. r*, but inside it can mean either r 
or r„, since the expression is already first order. For Eq. (16) to be true for all 9 we require firstly 
that 

2 Gm(r) 
(f>' = 4irGp-2n 2 , or to zero order 6' « ^ , (17a, b) 

and secondly, from the P2 coefficient, we obtain after some rearrangement 

„ 6a 2rm' f a! a\ , . 

a" - — + - + - = , 18 

r z m \r r z J 

where we only need the zero-order approximation (17b) to eliminate (j), since a is already 
first-order. We can therefore determine a by first solving the usual stellar structure equations 
to obtain m(r),m'(r), and then integrating Eq. (18) with these functions, subject to boundary 
conditions that a and a' are finite at the centre. These boundary conditions determine a up to a 
multiplicative constant. For at r = 0, we have rm'/m = 3. Hence from Eq. (18) a ~ B + C/r 5 , 
and so the central boundary conditions require C = 0, leaving only an arbitrary multiplicative 
constant. Consequently Q is determined unambiguously, since the expression (15c) for Q is 
independent of a constant factor in a. For a polytrope of index n in the range - 4.95, we find 
that Q as given by Eqs. (15c) and (18) can be approximated by the interpolation formula 

/ „\ 2.215 

' - a \ n.0245n-0.096ri 2 -0.0084n 3 



Q „ ^(l-^J e U-™-U.U 9 bn-U.UU«4^ ±1 . 5%r . m . s .. (19) 

(v) For a star distorted by the gravitational field of a companion, and not rotating, the calculation 
is the same as in (iii) except that, by reference to Equn (3), in the potential (13) we must replace 
S7 2 /3 by —Gm2/d 3 inside the parentheses; outside the parentheses we put 17 = 0. Also, of course, 
the angle on which P2 depends is measured now from the line of centres, rather than from the axis 
of rotation. By analogy with Eq. (15), 

m 2 rf 1 p m 2 r\ Q 

ai ~-^T^Q ' q ~ d3 l-Q ' (2 ° a ' 6) 



with the same structure constant Q as before. 
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(vi) Schwarzschild (1958) obtained, by a slightly different route, a result which in our notation is 

m 2 rf (3a — rat' 



comp _ i — /■oil 

9 " d 3 ^a + ra'A^ ' 1 j 

Although superficially very different from Eq. (20b) with Q given by Eq. (15c), it is in fact the 
same by virtue of the fact that 

— mr 2 (3a — ra') = r 2 (5a + ra')m' , (22) 

as can be verified by using Eq. (18) to eliminate a" from the LHS of Eq. (22). Thus 

= ( ^-ra> \ = 1 h ha + ra r )r 2 dm (23) 

V 5a /r=n bm\r^a\ J 

(vii) We can obtain a crude but illustrative approximation to a by assuming that (a) the mass 
is concentrated at the centre, and (b) the quadrupole is concentrated at the surface, where the 
distortion is greatest. Then the potential inside the star (r* < r±) - cf. Eq. (12) for outside - is 
given by 

Gmi Gy ot r 2 P 2 (cos fl) 1 2 2 _ 
r r\ 3 r 3 

The condition that <\> is constant on any interior potential r* = const, implies that 



m „ ^± + ' ^ ' + ?nV ~ ^VP 2 (cos 9) . (24) 



,/ n Gmi 1 2 2 / Gmi Gq rot r 2 1 2 2 \ n , ^ 

-0(r*) « - + -ft 2 r* 2 + a + — d — -Sl 2 r 2 P 2 (cos6<) (25) 

is independent of and hence that 

a = r , 3 (J^_jr\ s MiL_j^ . (26) 

\3Gmi m ir \) \3Gm 1 m x r\) K ' 

Eq. (14) is just Eq. (26) evaluated at r* = r±. Putting a oc r 3 in Eq. (15c), we obtain an 
approximation for Q which we can evaluate numerically for polytropes and check against Eq. (19): 

.1^1 „ 3 ( 1 _n\ 2 - 205 e -0.437n+0.066n 2 -0.023n3 ± 3% r . m . g . 



_ .^^^ _ r M__ e -0.437n+U.0bbn"-0.l«3n- ± 3% r . m . s . . f 2 7) 

* 5 rfmi 5 V 5/ V ; 

The approximation atxr 3 that we have obtained can be seen to satisfy Eq. (18) in the limit of 
a centrally-condensed star, where m! = everywhere except at a central point. The full solution 
in this case is a = Br 3 + C/r 2 . We cannot assume that C is negligible despite the fact that it 
leads to singularity at the origin, since m' is also singular at the origin. However in practice, C 
does become negligible in polytropes as n — > 5. 

Note that the approximation a oc r 3 of Eq. (26) gives zero if we use the differential form of 
Eq. (23), as against Eq. (27) if we use the integral form of Eq. (23). The differential form requires 
us to obtain a more accurate solution for a than does the integral form, as is not unusual. 
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In the extreme that the star is of constant density, rm'/m = 3 throughout and Eq. (18) 
gives a = B + C/r 5 . This means that a is constant, since the solution must be non-singular at 
the origin. This is just the well-known case of 'liquid' stars. For such stars (Maclaurin ellipsoids, 
Chandrasekhar 1969) all the equipotentials are similar ellipsoids with eccentricity e where 

tfrt = !^EZ {(3 _ 2e 2 )sin -i e _3 ev ^ } „ | e 2 (2g) 
Grai 2e d ' 5 

The quadrupole moment of a uniform ellipsoid is 

,rot = _l mir 2 e 2 (1 _ e2) l/3 ^ . (29 ) 

This 'agrees' with Eqs. (15b) and (27) in the case that p =const. and hence Q = 3/5. The 
agreement is providential, however, since we ought to use Eq. (15c) with a =const., whereas Eq. 
(27) assumed a oc r 3 ; but both expressions give the same answer if p =const. This agreement, 
providential or not, coupled with the genuine agreement for centrally condensed stars (n — ► 5), 
appears to suggest that approximation (27) might be good enough over the whole range of models 
from uniform density to centrally condensed. However if, for polytropes, we compare the Q of Eq. 
(19), obtained from the direct numerical solution of Eqs. (18) and (15c), with Q obtained from 
Eq. (27), using the fitting formulae as given, we find that the disagreement can approach ~ 50% 
for the important range of n ~ 1.5-3. We therefore recommend against using approximations 
(26) and (27), convenient as they are for illustrative purposes. 

(viii) Putting together Eqs. (15) and (24), and recalling that the symmetry axes of the two 
quadrupoles are Q and d, we see that the quadrupole tensor of a star in a binary can be written 

QH = - - tfdij) + (Midj - d 2 5 lj ) , A = , (30a, b) 

where Q might be given either by Eqs. (15c) and (18), or by approximation (27). 

The quadrupolar distortion of Eq. (25) produces an acceleration, in addition to the gravity of 
two point masses, which we call fQD- This non-dissipative term is 

Gm 2 d i djq i i m 2 ^(Q.d) 2 m 2 AQ? Gm\A 
pf QB = -»= j5 = 2^^ + ^^ + ^^ < ( 31 ) 

and so the orbital equation is 

■ ; Gmd „ . . 

+ f Q D , (32) 



d 3 

where 



c m 2 A 

tQD = 

A* 



5(0.d) 2 d ft 2 d Q.dQ 6Gm 2 d 
2dP ~2dF d? dT~ 



(33) 
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3. Dissipation of the Equilibrium Tide 

We continue to suppose that *1 is extended while *2 is a point mass. In the frame which 
rotates with *1, the quadrupolar tide will in general be time-dependent: *1 will be continually 
changing its shape. If the star is not perfectly elastic we expect a loss of total mechanical energy, 
but no loss of total angular momentum for the system as a whole. 

In the previous § we had an acceleration fQD, Eq. (33), which is derivable from a potential 
(5$, Eq. (31). Under prevailing circumstances this conserves total energy, (including rotational 
energy), as we show shortly. But if there is in addition a slow dissipation of energy by tidal 
friction, we will have an extra acceleration fpF, say. Writing 

d = -^-I V <5d> + f TF , (34) 
d 6 n 

we can see that total angular momentum 

H = //dAd + itt (35) 
is conserved if the couple on *1 is given by 

ISl = dA(V»-pf TF ) • (36) 

The total energy, 

E = ^d.d-^^ + 5^(d,n.d) + lm 2 , (37) 

may change not only because d varies but also because Ji can vary, for instance as a result of 
precession. We find however that E = (in the absence of tidal friction) provided that 5$ 
depends only on d and Q. d: 

■ •• Gmd 1 
E = /id.{d + — =- + -V<f$} + fi.(d £$' + /«) 

d 6 fi 

= fid. f TF + j A (V<5$ - n f T F).(d 5& + IQ.) 

= n(d- n Ad). frF + V<5$.0 Ad (38) 

where we have used Eqs. (34) and (36) for d and Q. The quantity 5& means the partial derivative 
of 5& with respect to Q.d. For 5$ a function only of the quantities d and fi.d, we see that V<5$ is 
entirely in the plane of d and £1, and hence that the last term in Eq. (38) vanishes: 

E = fi— .fpF , since d = — + O A d , (39a, b) 

where d/dt is a derivative in the frame rotating with *1 and d is the derivative in an inertial 
frame. 
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When there is dissipation, probably the simplest assumption we can make is that the rate of 
loss of energy is some positive-definite function of the rate of change (as seen in the frame rotating 
with the star) of *l's shape, e.g. of its quadrupole tensor since this determines its shape to lowest 
order. We therefore write 

dqij dqij 



E 



-a ■ 



dt dt 

where qij is given by Eq. (30) and a is a dissipation constant intrinsic to *1 (dimensions 



(40) 



dd/dt, so that 



2 t 1 ). In the frame that rotates with *1, ft is approximately constant while d varies at rate 



dq 



dt 



v _ 



d m<iA 
dt 2d 5 

3777,2^4 
2d 5 



(Sdidj — d 2 5ij) 



7 ddj 

di^ + d 
dt 



■'j 

dd, 
°~di 



+ d- Slj 



5dd , , 

dm d ^ 



Squaring this, after some manipulation we obtain 



9am\A 2 



E 



MY 
dt) 



+ 



dd\' 
~dt) 



9amlA 2 dd 
2d 10 ~dt 



2dd. 



dd 

~dt 



,dd 



(41) 
(42) 

(43) 



Comparing this with Eq. (39a) we see that a consistent expression for the acceleration fpF due to 
tidal friction is 

9am 2 . A 2 r . . dd „dd~ 



fpF = — " 



2fid w 



2ddr ft +d2 Tt 



9am\A 2 
2^d w 



3dd.d+ (h- Ud 2 ) Ad 



(44) 
(45) 



The mathematical form of this expression for the perturbing force due to tidal friction is 
exactly the same as that obtained by the more traditional model in which it is assumed that the 
tide, which would be directed along the line of centres if it were able to adjust instantaneously to 
the potential, is offset by some small lag time that is proportional to the timescale of dissipation. 
It is not clear that this lag time should be independent of orbital phase, but such an assumption 
is justified on the basis of our own model. 



4. The Effect of a Perturbing Force on the Orbit 

Using the relative position vector d = di — d2, and mass m = m\ + m2, the relative motion 
of a binary subject to (i) Newtonian point-mass gravity, and (ii) an additional acceleration f (the 
perturbing force per unit reduced mass \i = mim-z/m) is given by 

•■ Gmd „ , . 

d = ~^- + f. (46) 

A remarkably straightforward way to analyse the effect of a general (small) force f on the otherwise 
Keplerian orbital solution is by way of the LRL vector (Heggie, personal communication, 1996). 
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This has many applications apart from the present one (Eggleton 2000, in preparation), and so we 
describe it in a slightly more general way than is strictly necessary for present purposes. 

Define £ (Keplerian energy, per unit reduced mass), h (angular momentum, similarly) and e 
(LRL vector) by 

1 • • Gm , , • Gmd . t _ . 

£ = -d.d , h = d A d , Gme = d A h — . (47a - c) 

2d a 

Note that £ is not the total energy, if f ^ 0, but only the part that is kinetic plus Newtonian 
point-mass energy. We can see, after some manipulation in the case of e, that 

£ = d.f , h = d A f , Gme = f Ah + dA (d Af) . (48a - c) 

Hence £ , h and e are all constants of the motion if f=0. Using a standard parametrisation of the 
Keplerian orbit, for example either of the two parametrisations (51) - (55) below, we find that the 
vector e defined in Eq. (47c) is in the direction of periastron, and has magnitude equal to the 
eccentricity (thus justifying belatedly the choice of name e). In this section we show how £, h and 
e can be estimated by averaging the LHS of Eqs (48a) - (48c) over a Keplerian orbit, supposing f 
is small. 

Even if f is not small, auxiliary variables a (semimajor axis), b (semiminor axis), I (semilatus 
rectum), uj (mean angular velocity) and p (period) can be defined in terms of £,h, e in the usual 
way: 

a = — , b = ay/1 — e z , / = all - e , uj = — = — . (49a — a) 

2£ ab p 

For general f, and not just f « 0, four standard relations can be shown to be satisfied identically: 

2h 2 £ + G 2 m 2 {l-e 2 ) = , e.h = , h 2 = Gml , uj 2 = Gm/a 3 

(50a - d) 

Thus even when h, l,oo and a are continuously changing because f / 0, the orbit can be perceived 
as always 'instantaneously Keplerian'. For example, the instaneous period and semimajor axis 
always satisfy Kepler's third law, Eq. (50d). 

If f is small, we can estimate its effect on £,h, e by averaging over time the RHS's of Eqs 
(48) in a Keplerian orbit. This is done most easily by writing the Keplerian orbit in Cartesians 
(origin at focus, e in the 1-direction, q = h A e in the 2-direction, h in the 3-direction) using one 
or other of the following parametric forms: 

d = 77 (cos 6, sin 6, 0) = (a cos 4> — ae, b sin (f>, 0) (51) 

1 + e cos 6 

d = (— sin#, cos 6* + e, 0) = — (— a sin <p, 6 cos (f>, 0) (52) 

I 1 — e cos <p 

I 2 d9 

wdt = — j— --k = (1 -e cos (p)d(f) (53) 

ab (1 + ecosfc*)^ 



-12- 



l 



d 



(1 + ecosfl) 
tuabe 



a(l — ecos</>) 



sin 6 



uae sin ( 



(54) 
(55) 



I 1 — e cos < 

Various scalar, vector and tensor functions of d can be averaged over time using these 
parametrisations, 6 being more useful if the function contains a substantial negative ( < — 2) 
power of d and <fi otherwise ( > — 1). Since e, q = h A e and h are not unit vectors, we use 
bars to indicate the corresponding unit vectors; then for example di = dei cos 8 + dljiSinO. Some 
specific averages, expressed in terms of polynomials I n j(e), n, I > 0, defined below, are as follows: 



1 1 fP dt 

< d n + 2 > ~ n L ^+2 



1 fP dt 1 f 27r 

p Jo d n + 2 ' ' abl™ Jo 



d6 



(1 + ecos(9) T 



abl T - 



I 



?1,0 



^ jn—\ t 

< a > — a l n n 



n—1 



didi to 2 ab 



< d^ 2> = 7 ^{(l + e 2 Kn,o + 2e/ n , 1 } 
< d n didi > = uo 2 a n+2 {I nfl - e/ n> i} 
io 2 abe 2 



d 

< ~i — To > 



< 



d n + 2 

di 

d n+3 



ln+2 



{ln,0 — In,l} 



> = 



1 

P 



p cos 6 dt 
o H^ 2 



+ Qi 



p sin 0dt 



o d n + 2 



abl n 



< d n l di > = -a n {I rhl + el nfi }ei 



< 



d n+2 



> 



dd. 



toe 



{el nfi +In,l}qi 



d n+3 
didj_ 



didjdk 
d n + 5 

didjdk 



I 
1 
P 



(56) 

(57) 

(58) 
(59) 

(60) 

(61) 
(62) 
(63) 

(64) 



p cos 2 6 dt 



d n+2 
1 



p cos sin dt 



o 



-{-fn,2eiej + (J nj0 - In^lilj] 



p sin 2 tft 



o 



- In,3)(eiqjq k + + ^e fc ) + /n^e,-^} 



< > = p+r{( 7 n,3 - /n,i)(ci?je fc + g^j-efc) + (/ n ,i + el nfi )q iqj q k 

+(I n ,3 + el ny 2){e{ejq k - mjlk)} 
The polynomials /„, m , all positive, are defined by 

r 27T 

2vr 



m (e)= {** (1 + e cos 0)" cos m • 
Jo 



(65) 
(66) 

(67) 
(68) 
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They are easily evaluated from 

^ ~ dip (2m)! 



f 2 dtp (2m) ! 

h,2rn = j COS m 1p = 2 2m( m ^2 ' ^0,2m+l = , i"n+l,m = In,m + el n ,m+l 

(69) 

Clearly 

Jn,m(-e) = (-l) m /„, m (e) . (70) 



2vr 

In this paper the results that we need all come from the 0-parametrisation, but we include the 
0-parametrisation in Eqs (51) - (55) for completeness. 



5. The Effect of fgrj and f T p on the Orbit 

We can now use some of Eqs. (56) - (67) to average the RHSs of Eqs. (40a) - (40c) for the 
particular forces fQD, frF of Eqs. (33) and (45). We firstly consider the simpler case where the 
spin is parallel to the orbit, and then the more general case. 



5.1. Stellar Spin Parallel to Orbit 



We first specialise to the case Q. || h, which means that Q.d = 0. This puts to zero the most 
complicated term which appears in Eq. (33). Averaging £ by Eq. (48a), and consequentially a,u 
by Eqs. (49a), (50d), the potential force fQD gives zero. The term £tf of Eq. (45), using the 
averages (56) and (60) with Eqs. (68), (69), gives 



2Co 
3 to 



a 
a 



£ 
£ 



< d. f TF > 



9am 2 A 2 
2fj£ 



9am\A 2 u 2 ab 
2^£ W~ 



d 8 

2 m 2 

3e (^6,0 - h,2) + h,o — — r h,o 



3 < > +h 2 < 4tt > -fth < -4 > 



d 10 



d 8 



1 



31„2 , 255^4 , 185 .,6 , 25„8 
c -r 16 e -r 64 e 



1 + f e 2 + 



Q 1 + f e 2 + f e 4 + ^e 6 



[I 



2)15/2 



(1 



=2\6 



where the timescale £tf is defined, using Eqs (25), (49) and (50), by 



ixa 
9am 2 , A 2 



1 a m\ 



9(J rj u mim \ Q 



(71) 



(72) 



Similarly we average h with Eq. (48b). Once again fQD does not contribute, this time because 
it is purely radial if Q.d = 0; but fr F , using Eqs. (56), (68) and (69), gives the orbital average 



< d A f TF > 



9am\A 2 , 

7T- h 

2/x 



i n i 

< d 8 > ~h < d® > 
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9am%A 2 h 
2fiabl 6 



'6,0 



M 2 
ujab 



i.e. 



h 
h 



2£tf 



1 + f e 2 + f e 4 + ^e 6 Q 1 + 3e 2 + |e 4 



!)13/2 



(1 



,2\5 



(73) 
(74) 



We can also average e using Eq. (48c). We consider separately the terms from fq>D and fpF- 
Since in the present case (J) || h, Sl.d = 0) fQD is parallel to d, we obtain from Eq. (33) the orbital 
average 



Gme QD = < f QD > Ah = ^^hA < ^ > +— -^hA < ^ > 



(75) 



and so using Eqs (49), (61), (68) and (69), 



6qd = Z h A e 



Z = 



mA 
2mia 5 iO 



Q 2 , 30 m 2 (j 2 1 + |e 2 + \e A 
m (1 — e 2 ) 5 



(1-e 2 ) 2 



+ 



(76) 



The vector e is in the direction of the line of apses, and thus Z is the rate of apsidal motion, i.e. 
the rate at which the line of apses rotates about the direction of h. 

The effect of fpF on the average of e is to produce 5 terms: 

dd , , 9 d d d , d „ /mN 

e TF oc 3 < -p > Ah + h 2 < > -Qh < ^ > + < ^ > Ah - < ^ > AO . (77) 

By reference to Eqs. (64), (61) and (63), all these terms are ultimately parallel to e, i.e. they 
circularise the orbit without (further) apsidal motion. In this case we can obtain e more simply 
than e by differentiating Eq. (50a), and using Eqs. (71) and (74): 



e 
e 



1-e 2 



'£_ /T 
2E + h, 



9 



1 + 15 2 , 15 4 , 5_ 6 
(!_ e 2)13/2 



lin 1 + f e 2 + |e 4 
18u (1 - e 2 ) 5 



(78) 



We can obtain the rate of change of the intrinsic spin S7, using the constancy of H in Eq. (35): 

Cl fih h 7 

Q ~ ~7tt h ~ 2t T F 

where h/h comes directly from Eq. (74). The factor 7, the ratio of orbital to spin angular 
momentum, is usually large. As a result, the ratio Q/iv tends rapidly to a 'pseudo-synchronous' 
value (Hut 1981) which depends only on e and is obtained by equating to zero the contents of the 
brackets in Eq. (79). Then, on a slower time scale, e — ► and — ► u. 



1 + f e 2 + f e 4 



T 16 K 



n 1 + 3e 2 + |e 4 



(1 



g2)13/2 



10 



(1 



3 2\5 



7 = 



(79) 
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5.2. The Non- Aligned Case 

For the general case where there are rather more terms, but the problem is quite 

tractable. Eqs. (48) give £ , e and h in terms of the perturbing force; we average these over a 
Keplerian orbit, using the forces £qd and itf of Eqs. (33), (45). The recipes of Eqs. (56) - (70) 
allow us to express the results in the form 

8/8 = U (80) 

e/e = -Ve + Zq-Yh (81) 

h/h = Ye- Xq-Wh , (82) 

where the coefficients U, ...,Z may be somewhat complicated as functions of e, via the integrals 
In,m{e), but are otherwise simple functions of a, b, I, u>, a, .... The coefficients (X, Y, Z) are in fact a 
vector, the angular velocity of the e, q, h frame relative to an inertial frame: in § 5.1, X = Y = 0, 
and Z gives the term for apsidal motion in Eq. (76). The evaluation of U, Z is in our experience 
most easily done by splitting £1 into a component Oy = ii.hh parallel to h, and the corresponding 
perpendicular component £l±. Then we split each of the forces f*QD and fpF into two parts. One 
part, say t)f*QD or <5itf, is defined as 



_ m 2 A 
oiqd = 



5(o ± .d) 2 n ± .dn 

2d7 dF~ 



_ 9am 2 2 A 2 OjAd 

= —2H ¥~ • (83) 



These vanish when Q± = 0. The other part is the remainder. This remainder gives the same 
values for U = £/£, V = e/e, and W = h/h as Eqs. (71), (74), (78), except that f2 is to be 
replaced by fiy. It also gives the same apsidal motion Z as Eq. (76), without replacing O by Q.\\. 
The extra terms (5fqp, ^tf) brought into h by non-alignment are all in the plane of e, q, and 
contribute only to rotation. 

It is helpful to use Euler angles rj,X,i/j say to determine the orientations of e,q, h relative to 
an inertial frame, say E, Q, H, a suitable choice for H being the total angular momentum vector 
ofEq. (35): 

H = fih + m . (84) 

E and Q are arbitrary, provided they make a right-handed orthogonal set with H. The 
transformation from E, Q, H to e, q, h is the product of three successive simple rotations: by X 
about H, by r\ about (new) E, and by ip about (newer still) H, which now coincides with h. This 
gives 

e = (cos X cost/ 1 — sin X sin cos r/, sin X cos + cos X sin ^ cos 77, sin 77 sin (85) 
q = (— cos X simp — sinX cos ip cos 77, — sinX simp + cos X cosip cos rj, sin 77 cos ip) (86) 
h = (sin 77 sin X, — sin 77 cos X, cost/) , (87) 
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where we take the 1,2,3 - axes in the directions of E, Q,H. By differentiating Eqs. (85) - (87), it 
is straightforward to show that the angular velocity (X,Y,Z) of Eqs. (81), (82) relates to r),X,ip 
by 

X = fj cos tp + X sin tp sin 77 



Y = — 77 sin ^ + Xcos ^ sin 77 (89) 
Z = ip + X cos 7] . (90) 
From Eqs. (84) - (87), the components of in the directions of e, q and h are given by 

Q e = J2.e = y (H — /_th).e = Slo sin 77 sin -0 = £l±simp , (91) 

Q q = J2.q = Qq sin 77 cos ip = Q±cosi[) , (92) 

$7 = O.h = S7o cos 77 — fih/I (93) 

where £Iq, Oj_ have values 

IQ = H , ± = ^ sin 77 . (94) 

To obtain h and e from Eqs. (48b, c) and (45), we use most of the averages (56) - (67); 
we have to average quantities which range from scalars to tensors of rank 3. After substantial 
manipulation of the orbit- averaged expressions, which includes a pleasantly surprising amount of 
simplification, we obtain 

, . . Q± l + 3e 2 + fe 4 + fe 2 (l + ie 2 )cos2V> 
77 = Xcos^-Ysm^ = -— 8 {1 lJ )5 6 > 05) 

X sin 7/; + y cos t/> _ Am 2 ^o^l\\ 3fi e 2 (l + \e 2 ) sin 2ip 

X ~ Ih^i ~ ~2/W(l -e 2 ) 2 ~ 8^ (1 - e 2 ) 5 ' (%) 

and 

. Am 2 {n\-\^l) 15G m 2 2 A 1 + |e 2 + ±e 4 

t/^ + Xcost? = Z = 1 2 + § -> -f— 97 

2/icja°(l — e^J z /jwgt (1 — e z )° 

We see that in the absence of tidal friction (a = = 1/*tf) we have steady precession, with 77 = 
const, and 

^4777,2^0^11 , . 

X = — r- ^777 = const. (98) 

2/Jwa 5 (l - e 2 ) 2 1 ; 

We also see that Eq. (97) for apsidal motion is remarkably little different from Eq. (76) derived 
for the case of spin parallel to orbit. We emphasise that Eqs. (95) - (97) do not involve any 
approximation about the smallness of the inclination 77, or the eccentricity e. 

Our complete set of equations for h, e, 77, X and ip are then Eqs. (74) and (78) - with £1 
replaced by Qu - and (95) to (97). The ancillary variables a,(jj,£lu,Q± and Qq are given in terms 
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of h,e,7),H by Eqs. (49), (50) and (93), (94). A and t TF are given by Eqs. (25) and (72). We 
suppose that mi, m2, r±, Q, a, H, I are all given constants. 

The quantity defined in Eq. (94) is usually rather large: the angular velocity that the star 
would have if all the orbital angular momentum were converted into spin. Consequently rj, the 
angle between the orbital angular momentum and the total angular momentum, is usually rather 
small. The inclination i of the rotation axis of the star to the orbital axis is given by 

U± = Qsini = Oosin?? . (99) 

6. The Tidal Velocity Field and the Dissipation Constant 

In the above derivation the dissipation constant a was assumed to be given. However it is 
presumably determined by the internal physics of *1. The most important ingredients are likely 
to be (a) viscosity, with turbulent viscosity expected to be far more significant than molecular 
viscosity, at least in the regions that are turbulent, and (b) thermal diffusion. We shall assume 
here that turbulent viscosity is always the dominant dissipative agent, even in stars which are of 
early spectral type and normally assumed to have non-convective envelopes. Partly this is because 
even hot stars have at least one and usually two surface convection zones, driven by helium 
ionisation and/or bumps in the opacity. These zones may be small in mass but have considerable 
convective velocities, which presumably spill over into adjacent regions; they are also positioned 
close to the surface where they can contribute best to tidal dissipation. 

A further reason for expecting turbulence in the surface layers of upper MS stars is that they 
are generally rapidly rotating, and the rotationally-driven circulation currents, even though small 
in the interior, can become significant in the layers near the photosphere. In fact, a somewhat 
literal interpretation of rotationally-driven circulation (§ 7 below) gives a velocity field which is 
singular wherever the entropy gradient is zero (i.e. at the upper and lower edge of any convection 
zone) as well as at the surface where p = 0. Such singularities are unphysical, and are presumably 
rendered finite by the turbulent viscosity that is driven by the large shear in these layers. We 
expect that the upshot will be an outer envelope that contains significant turbulent viscosity, and 
consequently it seems reasonable to assume that the main agent of dissipation is this turbulent 
viscosity. Molecular viscosity is normally negligible in comparison, and so too is thermal diffusion. 
We return to this discussion in § 8. 

Note that in early-type stars, turbulent friction in the convective core may contribute to 
turbulent viscosity. This contribution is not expected to be very efficient (Zahn 1966). For recent 
discussions, see Zahn (1989), Goldman & Mazeh (1994), and Goodman (1997). 

In addition to the circulation velocity field (§ 7) there is, in the case of eccentric orbits and/or 
non-corotating stars, a time-dependent tidal velocity field driven by the time-dependent character 
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of the distortion. This field can be determined by using (a) the continuity equation 

^ + V.pv = , (100) 
and (b) the constancy of p on equipotential surfaces, i.e. the fact that 

9 = P( r *) > r * = r + raP2(cos9) . (101a, b) 

Working in the frame which rotates with *1, we need only be concerned with the companion- 
induced distortion, so that a±, the surface value of the distortion parameter a, is given by Eq. 
(20a). We are entitled to ignore a second contribution, from rotation, to the asphericity of in 
Eq. (101b), because our analysis below depends only on dr*/dt, and the rotational contribution is 
of course (nearly) constant in the frame that rotates with *1. 

Let 

F = r 2 P 2 = ^(k.rf-K 2 , so that V 2 F = , r.VF = 2F , (102) 

where k = d/d. Since d is time- varying, both a and k depend on t, the former because a oc 
(Equn 20a). Then with p as a function of r* only, and viewed as a function of r,t, we obtain 

dp_ _ dp_dr\ _ dp fdaF | 3aG \ _ 3a dp / 1 ^ | 

9t dr* 3t dr* \dt r r J r dr* \ d dt / 

where 

G = -— = k.r— .r , so that V 2 G = , r.VG = 2G . (104) 

3 at dt y ' 

F and G are obviously orthogonal harmonic functions of degree 2. Consider the velocity field 
given by 

Then Eq. (100) is satisfied to first order, provided that 

dp(3 a dp . 1 f r dp 

i.e. p = / a—dr . (106) 



dr ct\ dr pct\ J ri dr 

The lower limit in the integral comes from the boundary condition that the outer surface (p = 0) 
is a surface which moves with the fluid, so that the velocity must be finite there despite the 
vanishing density. The function (5 is determined unambiguously by the structure of the star, via 
Eq. (18) determining a(r), and is well-behaved for polytropic (0 < n < 5) surfaces as p — > 0, 
despite the apparent singularity there. For the special case n = 0, i.e. uniform density, we have 
(3 = a/a>i = 1. Fig. 1 illustrates the behaviour of p, a and [3 as functions of r in four cases: 
polytropes of index n = 1.5 and 3, and ZAMS stars of masses 1 and 8 Mq. 
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In suffices, 



where = (Skikj 



ddt ' h lh~ ^t kj ' (107) 

The tensor Sij is symmetric and traceless. This allows us to calculate the rate of dissipation of 
energy due to the action of turbulent viscosity on this velocity field, and this dissipation in turn 
determines the strength of tidal friction. 



The rate-of-strain tensor is now seen to be 



dvi dvj 3ai / j3' x \ 



(108) 



We square this and average it over an equipotential (which at this level of approximation can be 
taken to be spherical), using the results 



4vr J 



X 'i X j d — ^i>j 



We obtain 



Now, 



1 

47T 



If r 4 

— / XiXjX k xidQ = — {SijSki + S ik 5ji + SuSjk) 



tUn = 9a 2 s 2 



6 



ldd\' 
ddt) 



+ 2 



/dkV 
\~dt) 



2 

d 2 



dd\ 2 

m) + 



\~dt) 



(109) 
(110) 

(111) 



and so the rate of dissipation of mechanical energy is 

1 



E 



- j pwlt^dV 



9a\ 



m\ 2 

~dt) 



7 



mi / 



dm 



(112) 



The parameters w, I are the mean velocity and mean free path of turbulent eddies. As we 
emphasised at the start of this Section, w might be due to shear in the rotationally-driven 
circulation velocity field, as well as to convection (of which there is normally some even in stars 
with supposedly radiative envelopes). Fig. 1 shows the how the weight factor in parentheses in 
Eq. (112), which we call 7(7"), varies throughout some polytropic and ZAMS model stars. 

It is gratifying (but not coincidental) that Eqs. (Ill), (112) contain the same functional 
dependence on d(t) as Eq. (43), despite being derived by a very different route. Equating (43) 
with (112) allows us to determine a as 



a 



wlj(r)dm , 7 (r) = [3 2 + \r(30 + ^r 2 f3' 2 



m\r\Q 2 

This in turn, via Eq. (72), gives an estimate for the tidal friction timescale: 



18mm2 



J wh{r)dm 



(113) 



(114) 
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Fig. 1. — The variation of 0.1 \og(p/p c ) + 1 (plusses), ajct\ (asterisks) and (3 (circles) with radius 
(arbitrary units); p c is the central density. Also shown are n/n + 1 (continuous line), where n is the 
local effective polytropic index, and finally (crosses) 0.1 log 7 + 0.5, where 7 is the weight factor in 
Eq. (113), i.e. the contribution which each layer would make to the dissipation rate if the viscosity 
were uniform, (a) 8Mq (b) IMq (c) n = 3 polytrope (d) n = 1.5 polytrope. Modest fluctuations in 
n/n + 1 near the centre in the two stellar models are due to the fact that this quantity was obtained 
by the ratio of differences in a tabular model, which had converged to a maximum-modulus error 
of ~ 10~ 6 . 
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We defer to § 8 a brief discussion of the possible significance of this estimate, but we note here that 
it does not depend on any of the following approximations and assumptions: (i) V A v = 0, or 
V.v = 0, (ii) a oc r 3 , or (3 oc r 4 , (iii) ft \\ cj, or e ~ 0. It depends only on the approximation that 
the quadrupole tide dominates, that it is small enough to validate a linearisation of the continuity 
equation (100), and that there can be assumed to be an isotropic viscosity. 



7. Rotationally Driven Circulation 

In the surface layers of rapidly rotating stars, rotationally-driven circulation could well be the 
second most important source of turbulence, after convection. For a homogeneous (i.e. zero-age) 
star, T as well as p, p is constant on equipotentials, and hence so is e, the rate of release of nuclear 
energy. Consequently we can define a nuclear luminosity L interior to an equipotential by 



dL 
dV 



pe . (115) 



Let us define 



K(<f>) = AirGm - 2fl 2 V = J V0.dE = J |V0|d£ (116) 
and note, from Eqs (lb) and (5), that 



V^f . (117) 



Using Eq. (6) for the volume element, we see that hydrostatic equilibrium can be written as 

_ldp_ _# d^dV_ _ Airr* 2 _ ( Gm 2Q 2 rA (4?rr, 2 ) 2 
~p~d^ ~ ~dV* dV~dV* " Td£/W<t>\ ' Vn 2 ~J /|V0|d£/d£/|V0| ' 



(118) 



The factor in parentheses simply cancels the first factor in the denominator to its right, by (116); 
we write it this way to show that the ratio on the right of the parentheses clearly differs from 
unity in second order if (p differs from spherical in first order, so that we can write 

1 dp _ d<j) _ Gm 20 2 r* 
p dr* dr* r* 2 3 

This is the first integral of Eq. (17a), keeping the first-order term which was neglected in Eq. 
(17b). 

Now consider the energy flux F, which in general is a combination of radiative and convective 
flux. In spherical symmetry, we usually write this as 
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where w is the mean velocity of convection and T5S is the mean heat excess of an upward-rising 
eddy. The mixing-length approximations for w and T5s are 



o 

W ~ 



TSS 



Tl 



dS_ 

dr 



dS 



(121) 



d\ogp\ 

where I is the mixing length, normally estimated by I ~ — dr/dlogp, and where the square 
brackets have the meaning [X] = max(X, 0). In a non-spherical situation the generalisation of 
the radiative term in the energy flux is obvious; and of the many possible generalisations for the 
convective term we choose 



4acT 3 



VT + p 



i.e. 



F = X(<£)V0 



3np 

X(<f>) EE - 



dS 



dlogp 

AacT 3 dT 
3np deft 



3/2 j 

dcj) V 



P 



T- 



dS 
dlogp 



3 ' 2 dr* 



(122) 
(123) 



Thus the equation of energy production and transport, taking into account the meridional 
circulation current v, is 

= pe-pTv.VS , (124) 
where v also satisfies the steady continuity equation 



V.pv = . 



(125) 



Note that we are considering only the rotational component to the distortion of the star, which is 
approximately constant over many orbits, whereas the companion-induced distortion fluctuates on 
the timescale of one orbit, and produces the tidal velocity field of the previous Section. 

We first establish that the circulation term carries no net energy across an equipotential 
surface, i.e. that / pTv.VS dV = 0, where the integral is over the interior of an equipotential 
surface. From thermodynamics and hydrostatic equilibrium, 



TdS = dU + pd- = d(U + -)--dp =d(U +- + <£) . 

P P P P 

Hence, using (125), 

J pTv.VS dV = J pv.V(U + - + (f>) dV = J (U + - + <f>)pv.dS . 



(126) 



(127) 



Since the expression in parentheses is constant on an equipotential it can come outside the last 
integral, and since / pv.dH = by Eq. (125) we have the result. Hence 



J F.dS = J XV(f>.dY, = J pedV = L . 
Since X is constant on equipotentials, we can write this, using Eq. (116), as 

L = X J Vcj).dY, = XK , i.e. X = L/K . 



(128) 



(129) 
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Using Eqs (6), (115), (117) and (129), Eq. (124) becomes 

L dK .„ ,,, d L dL ^dS „ , . 

KW + ™l4K = Hv-" T T^f • < 130 > 

and so, after some manipulation, 

where v± is the component of v in the direction of V(j>. The tangential component of v then comes 
from continuity, Eq. (125). 

The analysis of § 2, in particular Eq. (9), allows us to estimate the angular-dependent term 
in Eq. (131) for the circulation velocity driven by the uniform rotation, in terms of the distortion 
parameter a and its surface value (Equn 15a). It is 

— — / V4> dX - V4> / — — « - — * — P 2 cos (9 . 132 

|V0| y 7 \V(f>\ 3Gmi(l-Q)ai dr v ; v y 

In addition to the rotationally-driven circulation of Equn (132), there is in principle a circulation 
due to the part of the potential that comes from *2's gravity. But unless and until *1 is brought 
into synchronism, and the orbit circularised, this contribution will fluctuate about zero with the 
period of *l's relative rotation. It will therefore be insignificant until synchronism is reached. 

It is evident that the velocity field v± given by Eq. (131) is singular on surfaces dS/d(j) = 
and p = 0. For mathematical consistency we require that there be a current sheet of zero thickness 
at each of these singular surfaces, to return the flow, which diverges out of the polar point on these 
surfaces from both sides (or else converges into it, also from both sides). Presumably the flow in 
and near this sheet is turbulent because of the high shear, and this turbulence provides viscosity 
which renders the flow finite in practice. 

There are many complexities in the study of meridional circulation (Mestel 1965), including 
for example the influence of composition gradients, and the redistribution of angular velocity from 
the (supposed) uniform state. Nevertheless in the surface layers of rapidly rotating stars it could 
well be the second most important source of turbulence, after convection. For recent work on 
meridional circulation, cf. the papers by Zahn (1992), Urpin et al. (1996), and Talon et al. (1997). 



8. Conclusions 

We have set out a formalism for developing the equations governing the spin and orbital 
evolution of a binary subject to the 'equilibrium tide' model of tidal friction, along with the 
non-dissipative effect of the quadrupole moment of the tide. We believe it may be helpful to have 
all the analysis in one place, and the analysis of the orbit by the use of the LRL vector seems 
considerably more efficient, as well as more transparent, than by other means. It leads directly 
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to formulae applicable to any degree of inclination of the star to the orbit, and not just to small 
inclinations (Hut 1981). We also believe it may be helpful to have shown that the normal equation 
for tidal friction, based on the picture that the tidal bulge lags the line of centres by some small 
constant time, follows quite directly from a more physical model in which the dissipation is related 
to a positive-semidefinite function of the rate of change of the tidal deformation (as measured by 
the quadrupole tensor) in the frame that rotates with the star. The effective lag time r is related 
to the rate of dissipation a in our model by 

We have included a tentative determination of a (Eq. 113) in terms of a hypothetical 
isotropic viscosity due to the internal turbulent velocity distribution in the star. This is a much 
more uncertain area, since the internal dynamics, particularly but not exclusively as it relates to 
turbulence, is far from clear, and very different formulations have been presented in the literature 
(Zahn 1978, Scharlemann 1982, Campbell & Papaloizou 1983, Tassoul k Tassoul 1990, Rieutord 
& Zahn 1997). Regarding convectively driven turbulence, we note that our weight factor 7 in Eq. 
(113) is not especially small in the convective core of the 8Mq star, being ~ 0.003 — 0.01 (Fig. la). 
This is larger than is expected with the rather inaccurate approximations a oc r 3 , (3 oc r 4 , 7 oc r 8 . 
We propose in the future to evaluate a from convection in a variety of stellar models, but a 
provisional estimate suggests a timescale txF ~ 30 yr, for 77 ~ a, mi ~ 7112 ~ 8Mq. 

Some further turbulence may be due to rotationally driven circulation (§ 7), which, while 
presumably not actually singular at the surface and at the boundaries of convection zones (of 
which there are often one or two near the surface even in stars with 'radiative' envelopes), may 
be more significant there than is usually supposed. Another source might be differential rotation, 
itself brought about because tidal friction is likely to operate most strongly on the outer layers. 
Differential rotation leads to shear which should lead to turbulence. 

It has long been suspected (Maeder 1975) that stars have larger convective cores than can 
be accounted for by 'standard' assumptions about convective stability. Although several different 
models for this have been proposed, they are not consistent with each other (Eggleton 1983). 
But even in the absence of a definitive model, one can compare theoretical models containing 
an artificial amount of 'enhanced mixing' with those (relatively few) observed binaries, the £ 
Aur binaries, which have (a) a highly evolved (red or orange supergiant) component, and (b) 
fundamental data (masses etc.) accurate to better than 5%. Schroder et al. (1997) show that £ 
Aur stars give best agreement with theoretical models that mix cores of about 40% more mass 
than 'standard' models at ~ 8Mq. Thus there is almost certainly more motion taking place in 
stars than the simplest theoretical concepts allow, and this means that an estimate such as Eq 
(113) for a is bound to be a lower limit; by how much remains to be determined. 

Although there is much disagreement in detail about the kinds of velocity fields to be expected 
in stars, including disagreement about the basic assumption of hydrostatic equilibrium from which 
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all the calculation of §s 2-7 stem, we feel that at least this last assumption is reasonably secure at 
the level which we require. Velocities would have to be sonic throughout the star for hydrostatic 
equilibrium to be seriously modified. Perhaps they may be this high very close to the photosphere, 
in very rapidly rotating stars and/or stars with surface convection, but hardly in deeper layers. 
Our Eq. (113) gives quite substantial dissipation from the deeper layers, including the core. Our 
determination of a in terms of wl is tentative given that the distribution of turbulence is not 
well-known, but we believe our answer is 'definitive' otherwise, particularly since it leads to just 
the same functional form for dissipation as does the earlier, more general, treatment of § 3 which 
makes no reference to a specific agent of dissipation. 

We believe it may be better for the present to think of a as a parameter that might be 
amenable to observation. We consider briefly three situations where a, or at least limits to a, 
might be measurable. 

(i) PSR 0045-716: This radio pulsar in the SMC (51.2 d, e = 0.808; Kaspi et al. 1994) is unusual 
in that the star it orbits, an apparently normal B1V star, is not a Be star and therefore there 

is little or no accretion to perturb the orbit in an erratic fashion. The highly eccentric orbit is 
quite wide so that even at periastron the stars are not very close: the separation between the 
two stars at periastron is approximately 4 times the radius of the companion B-star (Kaspi et al. 
1996). The variation of the orbit with time is already measurable: the orbital period divided by 
the time derivative of the orbital period is 5 x 10 5 yr (Kaspi et al. 1996). This variation may be 
due predominantly to tidal friction. In a reasonable interval of time rates of change of period, 
eccentricity and inclination should all be measurable. Our model gives equations which make a 
definite prediction about the time- variation, apart from the dissipation rate a, and which may 
therefore be testable. This might not be quite definitive, since we need to know such things as both 
masses, and the orbital inclination, rather than just the (very well-determined) mass function; but 
in principle all the unknowns might be determined by fitting over a long stretch of observation. 
Lai et al. (1995) have made a strong case that the spin axis of the stars is inclined with respect to 
the total angular momentum of the system. The analysis of the dynamical tide for this system was 
carried out by Lai (1996), Kumar & Quataert (1997), and in more detail by Kumar & Quataert, 
(1998), and Lai (1998). In principle, all of e,r],X,ip might be measurable in time. 

(ii) A Tau: This is a well-known triple system in which the third body is very close to the 4d 
eclipsing pair (P ut ~ 33d; Struve & Ebbighausen 1956, Fekel & Tomkin 1982). The orbits are 
quite probably parallel, or conceivably anti-parallel (Soderhjelm 1975). In the latter case the total 
angular momentum would be close to zero. The third body should be injecting a fluctuating 
eccentricity into the close pair, on a timescale of days. For point masses, this fluctuation would be 
of order 0.007, if the outer orbit is circular. Such a fluctuation in the eccentricity of a semidetached 
pair might be expected to have a dramatic effect on the rate of mass transfer, since the pressure 
scale-height in the atmosphere, relative to the radius, is probably at least ten times smaller. The 
absence of such an effect suggests that tidal friction is quite efficient at removing eccentricity; 
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but too-strong tidal friction could mean that the inner orbit shrinks rapidly compared with a 
reasonable estimate of its age on the grounds of nuclear evolution. Thus we might expect to find 
both a lower and an upper limit to the dissipation constant a from this system. 

(iii) Algol: The well-known triple star Algol (Tomkin & Lambert 1978, Hill et al. 1971, Labeyrie 
et al. 1974, Popper 1980) has a long-period orbit (P ut = 1-86 y) which, from VLBI measurements, 
is inclined at ~ 100° to the orbit of the inner eclipsing pair (Lestrade et al. 1993). In such a 
highly-inclined system, if the stars are treated as three point masses the inner orbit should cycle in 
eccentricity between e ~ and e ~ 0.98, on a timescale of ~ ^out/-^™ ~ 500 y, while preserving 
both inner and outer periods. It presumably does not do so, and this might partly be because of 
tidal friction. However, the quadrupolar distortion of the stars even in the absence of dissipation 
should lead to apsidal motion in the close pair, which would be quite rapid compared to the 
apsidal motion driven by the third body. Since the large cycles in eccentricity depend rather 
critically on the rate of apsidal motion driven by the third star, the extra apsidal motion of the 
quadrupolar distortion may prevent these cycles from being large; but we might still expect larger 
fluctuations on account of the inclination than we would if the orbits were parallel. It is not clear 
that there is a lower limit to a from the (presumed) absence of large cycles in eccentricity. But at 
the same time tidal friction must not transfer angular momentum from the inner pair to the outer 
pair so fast that the inner pair would have difficulty in lasting for a nuclear timescale, since the 
system has in fact existed long enough for considerable nuclear evolution. Thus there may still 
be an upper limit. Other close triples such as DM Per (Hilditch et al. 1986), VV Ori (Duerbeck 
1975), HD109648 (Jha et al. 1997) and C 2 CrB (Gordon & Mulliss 1997), with outer periods of 
lOOd, 119d, 121d and 250d, may also give upper limits. 

The first of the above three situations gives a reasonable chance of either confirming or ruling 
out the equilibrium tide model; all the initial conditions as well as a will have to be treated as free 
parameters, which might be determined by fitting over a considerable stretch of time. An analysis 
of the two triple systems (A Tau, [3 Per) is currently being undertaken. 

LGK and PPE are grateful to NATO for financial support with grant CRG 941288. PPE is 
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